Draft version February 2, 2008 

Preprint typeset using I4TgX style emulateapj v. 7/15/03 



HIGH RESOLUTION X-RAY SPECTROSCOPY OF THE FE K COMPLEX IN IC 4329A 
Barry McKernan 1 - 2 and Tahir Yaqoob 2,3 

Draft version February 2, 2008 

ABSTRACT 

We report the detection of complex Fe K line emission from a Chandra High Energy Transmission Grating 
Spectrometer (HETGS ) observation of the Seyfert 1 galaxy IC 4329A. The line emission is double-peaked, 
one peak centered at ~ 6.3 keV, and the other at ~ 6.9 keV in the source rest frame. When modeled by 
Gaussians, the lower energy peak is resolved by the HEG at > 99% confidence, whilst the higher energy peak 
is resolved at only < 90% confidence. The best-fitting widths are ~ 21 ,000 km s" 1 and ~ 4000 km s" 1 FWHM 
for the ^6.3 keV and ^6.9 keV peaks respectively. If the peaks correspond to two distinct emission lines, 
then the peak energies are redshifted with respect to the expected line energies of Fe I Ka and Fe XXVI Lya by 
at least 650 kms" 1 and 950 km s" 1 respectively. Alternatively, the Fe K line profile may be due to a single line 
from a relativistic accretion disk. In that case the inclination angle of the disk is required to be 24+j degrees, 
the outer radius constrained to several tens of gravitational radii, and the radial line emissivity flatter than r" 01 . 
Another possibility is that both peaks are due to distinct lines but each one relativistic ally broadened by a 
disk. In that case the lower energy peak could correspond to emission from Fe in a low ionization state, and 
the high-energy peak to Fe XXVI Lya emission. Then, the inclination angle is even less, restricted to a few 
degrees. However, the radial emissivity law is allowed to be steeper (~ r" 2 5 ) and the outer radius does not 
have to be fine-tuned. Yet another scenario is that the lower energy peak originates in a disk but the higher 
energy peak originates in more distant matter. The disk inclination angle is then intermediate between the 
last two cases but the emissivity is again required to be flat. We cannot rule out Fe XXV He-like absorption 
modifying the observed line profile. However, the data, and inferred emission-line parameters, are insensitive to 
the presence of a Compton reflection continuum. Including Compton reflection does, however, allow a steeper 
radial emissivity law for the relativistic line. Future missions such as Astro-E2 will be able to break a lot of 
the degeneracy in the physically distinct models that can all account for the Chandra data. Since IC 4329A is 
one of the brightest Seyfert 1 galaxies it should be a good astrophysical laboratory for studying the ionization 
structure of accretion disks around supermassive black holes. 

Subject headings: galaxies: active - galaxies: individual (IC 4329A ) - galaxies: Seyfert - techniques: spec- 
troscopic - X-rays: line - emission: accretion - disks: galaxies 



1. INTRODUCTION 

Some of the Fe Ka fluorescent emission line in active galac- 
tic nuclei (AGNs) is believed to originate in a relativistic ac- 
cretion disk around a supermassive black hole (see e.g. re- 
views by Fabian et al. 2000; Reynolds & Nowak 2003). Early 
studies of type 1 AGN using X-ray observations with ASCA 
revealed a wide variety of Fe Ka line widths and shapes (e.g. 
Nandra et al. 1997; Lubiriski, P. & Zdziarski 2001; Perola 
et al. 2002; Yaqoob et al. 2002). These results have been 
confirmed with the advent of the Chandra and XMM-Newton 
X-ray observatories, but a more complex picture is emerging. 
New, high resolution observations are revealing that a narrow 
Fe K line (< 10,000 km s" 1 FWHM) component is common 
in type 1 AGN (Yaqoob et al. 2001 ; Kaspi et al. 2002; Pounds 
et al. 2001; Turner et al. 2002; Fang et al. 2002; Reeves et 
al. 2002; Yaqoob & Padmanabhan 2004), often superimposed 
on top of a broad (~ 10,000-100,000 km s" 1 FHWM) com- 
ponent. Traditionally, narrow Fe K lines are associated with 
matter at least several thousand gravitational radii from the 
black hole (one gravitational radius, r g , defined as GM/c 2 ). 
However, there is increasing evidence to suggest that narrow 
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Fe K lines can originate in, or have a significant contribution 
from, the putative accretion disk (Petrucci et al. 2002; Lee et 
al. 2002; Yaqoob et al. 2003a). Narrow Fe K lines originat- 
ing from a disk can be interpreted in terms of truncated disk 
emission (e.g. Done, Madejski, & Zycki 2000) or in terms of 
a flattened radial line emissivity where the intensity per unit 
area falls off less steeply than ~ r~ 2 (Yaqoob et al. 2003a). 
Alternatively, the narrow Fe K lines could be due to local- 
ized emission from the disk (Turner et al. 2002; Yaqoob et al. 
2003a). There are, of course, other possible interpretations of 
the origin of narrow Fe K lines (e.g. Sulentic et al. 1998; 
Elvis 2000; Hartnoll & Blackman 2000). 

Earlier observations with ASCA also identified Fe K-shell 
emission due to highly ionized (H-like) Fe in two type 1 AGN, 
PG 1 1 16+215 (Nandra et al. 1996) and Ton S 180 (Turner et 
al. 1998; Comastri et al. 1998). Observations with the Chan- 
dra high energy transmission grating spectrometer (HETGS) 
and the XMM-Newton CCD detectors are now revealing Fe K 
emission due to H-like and He-like Fe in other AGN (e.g. 
Reeves et al. 2002; Dewangan 2002; Page, Davis, & Salvi 
2003; Yaqoob et al. 2003a). As more examples of highly ion- 
ized Fe spectral features are detected in the spectra of type 1 
AGN, it is clear that we need a better understanding of the 
complex ionization physics of the accretion disk (e.g. Nayak- 
shin & Kallman 2001; Ballantyne, Ross & Fabian 2001; Bal- 
lantyne & Ross 2002). 

Here we report on the spectral results of a Chandra HETGS 
observation of the type 1 AGN IC 4329 A, in which we ob- 
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serve a complex, broad, double-peaked line profile. Prelimi- 
nary results of this observation were presented in Yaqoob & 
Padmanabhan (2004) and we now present detailed modeling 
of the Fe K line profile at the highest spectral resolution yet 
available (FWHM - 1860 km s" 1 at 6.4 keV). IC 4329A is a 
bright Seyfert 1.2 AGN at z = 0.016054 (Willmer et al. 1991) 
embedded in a nearly edge-on host galaxy. It is one of the 
brightest Seyfert galaxies known, and is sometimes brighter 
than 3C 273 in the 2-10 keV band. IC 4329A is luminous, 
with L (2-10keV) ~ 4 x 10 43 ergs s" 1 typically, where we 
have converted from the luminosity given by Reynolds (1997) 
using Hq = 71 km s _1 Mpc _1 (Spergel et al. 2003) and qo = 0. 

IC 4329A was observed previously with ASCA (Mushotzky 
etal. 1995; Cappi et al. 1996), GRO (Madej ski et al. 1995), 
BeppoSAX (Perola et al. 1999), ASCA simultaneously with 
RXTE (Done, Madejski & Zycki 2000) and XMM-Newton si- 
multaneously with BeppoSAX (Gondoin et al. 2001). AGaus- 
sian model fitted to the ASCA observations (Done, Madejski 
& Zycki 2000) revealed strong (EW= 180 ±50 keV), broad 
(FWHM=42,960± 11,015 km s" 1 ) Fe K emission peaking 
near ^6.4 keV. The XMM-Newton observations (Gondoin et 
al. 2001) revealed a narrow Fe Ka line (FWHM < 8785 km 
s , EW= 43 ± 1 eV) originating in matter with a low ioniza- 
tion state (the line peaked at E = 6.42+[}5jjkeV). 

The paper is organized as follows. In |5]we discuss the ob- 
servation and data. In |3]we discuss fitting different physical 
models of the complex spectrum in the Fe K band. In |4]we 
discuss the implications of our results and present our conclu- 
sions. 




Observed Energy (keV) 

FIG. 1.— IC 4329A HEG spectrum between 3.0 and 8.0 keV binned at 
0.1 6A compared with the best-fitting power-law model, (a) Photon spectrum, 
(b) Ratio of the data to the model. 



2. OBSERVATION AND DATA 

IC 4329A was observed with the Chandra HETGS (Markert 
et al. 1995) on 2001 August 26 from UT 06:14:07 for a du- 
ration of ~ 60 ks. The Chandra HETGS consists of two grat- 
ing assemblies, a high-energy grating (HEG) and a medium- 
energy grating (MEG). We used only the combined ±1 orders 
of the grating data. The mean Chandra total HEG and MEG 
count rates were 1.158 ± 0.0051 ct/s and 2.244 ± 0.0067 ct/s 
respectively. The source flux showed little variability over the 



entire duration of the campaign. For example, for HEG plus 
MEG combined first-order lightcurves binned at 128 s, the 
excess variance above the expectation for Poisson noise (e.g. 
see Turner et al. 1999) was (3.4 ± 18.2) x 10" 5 , consistent 
with zero. Only single spectra per instrument were therefore 
extracted. 

We made effective area files (ARFs or ancillary response 
files), photon spectra and counts spectra following the method 
of Yaqoob et al. (2003b). A net exposure time of 59,087 s 
was obtained (including a deadtime factor of 0.0162). The 
HEG bandpass is ~ 0.8-10 keV and the MEG bandpass is 
~ 0.5-10 keV but the effective area of both instruments falls 
off rapidly at either end of the bandpass. Background was 
not subtracted since it is negligible in the energy region of 
interest. The HEG has twice the spectral resolution (0.012A 
FWHM) of the MEG (0.023A FWHM) and has a larger ef- 
fective area than the MEG in the Fe K band so our analysis 
will focus on the HEG data, since we are interested in study- 
ing the Fe K region with the highest spectral resolution cur- 
rently available. The MEG spectrum of IC 4329A is discussed 
elsewhere (McKernan, Yaqoob, & Reynolds 2004, in prepa- 
ration). 

3. MODELING OF THE FE K COMPLEX IN IC 4329A 

We used XSPEC v. 11.2.0 for spectral fitting to the HEG 
spectrum in the 3-8 keV band. This choice for the lower end 
of the bandpass avoids the known complexities in the soft X- 
ray spectrum due to absorption and the higher end ensures 
that background will indeed be negligble since the instrument 
efficiency area vanishes rapidly at high energies. Galactic ab- 
sorption of 4.55 x 10 20 cm" 2 (Elvis et al. 1989) has a neg- 
ligible effect in this energy range and so was not included in 
our modeling. There has been a well-documented degradation 
in the quantum efficiency of Chandra ACIS due to molecular 
contamination. However, the effect of this degradation is most 
pronounced at energies below ~ 1 keV and is negligible over 
the energy range discussed here so we do not include ACIS 
degradation in our modeling. The C-statistic (Cash 1976) 
was used for minimization in all models and all statistical er- 
rors quoted are 90% confidence for one interesting parameter 
(AC=2.706), unless otherwise stated. 

First we investigated how the 3.0-8.0 keV HEG data com- 
pared with a simple model consisting of a single power law 
plus intrinsic absorption, with the spectrum binned at 0.16A. 
We included absorption in this initial fit, just to ensure that 
it can be neglected in the bandpass that we used. The best- 
fitting power-law index for this model was 1 .66^ 05 and the 
absorbing column density was zero with an upper limit of 
5.1 x 10 21 cm" 2 . Figure 1(a) shows the HEG photon spec- 
trum of IC 4329 A with the model superimposed. Figure 1(b) 
shows the ratio of the data to the model. The single power 
law is clearly a good overall fit to the hard X-ray spectrum 
but there is a broad double-peaked feature in the Fe K band. 
The peaks occur at ~ 6.2 keV and ~ 6.8 keV in the observed 
frame, and the lower energy feature appears to be broader than 
the high-energy one. To investigate the Fe-K emission com- 
plex in more detail, we used spectra binned at 0.005 A, which 
corresponds approximately to the lcr HEG resolution. Since 
no absorption is required in the 3-8 keV band, we used only a 
power law for the continuum. Several models of the complex 
Fe K emission are discussed below. 

3.1. Double Gaussian Model 
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We added two Gaussian emission-line components to the 
power-law continuum model in order to fit the ^6.2 keV and 
~ 6.8 keV emission peaks. We did not include a Compton re- 
flection continuum in the fit since we will show post-facto that 
the emission-line parameters are insensitive to the p rese nce or 
absence of a Compton reflection continuum (see 33. 6> . The 
best-fitting parameters for this model are given in Table 1. 
Note that in Table 1 we show the velocity shift from systemic 
when the lower and higher-energy peaks are identified with 
Fe I Ka emission (at 6.400 keV in the rest frame) and H- 
like Fe XXVI Lya emission (at 6.966 keV rest frame) respec- 
tively. Even when we take into account the systematic uncer- 
tainty of the HEG wavelength scale (~ ±420, ±460 km s" 1 at 
~ 6.2 keV and ~ 6.8 keV respectively) 4 , the low-energy fea- 
ture (if Fe I Ka emission) is redshifted by > 650 km s _I and 
the higher energy feature (if Fe XXVI Lya ) is redshifted by 
> 950 km s" 1 (both are 90% confidence lower limits). If the 
higher energy feature is due to Fe XXV forbidden or Fe XXV 
resonance emission, then the blueshift is ~ 12,200 km s" 1 or 
~ 9, 220 km s -1 respectively, relative to systemic velocity. 
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FIG. 2. — Confidence contours of line intensity versus rest-frame energy 
for (a) the lower energy Fe K peak and (b) the higher energy pea k from a 
double Gaussian model fit to the HEG spectrum of IC 4329A (see i|3.1l and 
Table 1). The two-parameter joint confidence levels correspond to 68%, 90%, 
and 99% confidence. The dashed line in (a) denotes the rest-frame energy for 
Fe I Ka emission (6.400 keV) and that in (b) corresponds to the rest-energy 
of Fe XXVI Lya emission (6.966 keV). The contours in (a) indicate that the 
lower energy feature is redshifted from the Fe I Ka rest-energy at ~ 90% 
confidence. From (b) it appears that the higher energy feature is redshifted 
from the rest-energy of the Fe XXVI Lya emission line. 

The low-energy Fe K feature is moderately broad 

(FWHM=20, 825!'3 7 y° km s" 1 ), as we saw from 
Fig. 1, whereas the higher energy feature is narrower 
(FWHM=3680!^ km s _1 ). We estimated the significance 
of each feature by comparing the values of the C-statistic 
obtained from fitting a power-law continuum only and from 
adding the feature in question to the power-law only fit. The 
resulting values of AC are shown in Table 1. The addition 

4 http://space.mit.edu/CXC/calib/hetgcal.html 



of each Gaussian component adds three free parameters to 
a fit, so the lower energy Fe K feature was detected with a 
confidence level > 4 a and the higher energy feature was 
detected at only > 95% confidence. 
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FIG. 3. — Confidence contours of equivalent width (EW) versus FWHM 
for (a) the lower energy Fe K emission peak and (b) the higher energy peak 
from a double-Gaussian model fit to the HEG spectrum of IC 4329 A (see 
i!3. H and Table 1). The two-parameter joint confidence levels correspond to 
68% (innermost), 90% and 99% (outermost). 

Fig. 2(a) shows a two-parameter joint confidence contour 
of line intensity versus rest-frame energy for the lower en- 
ergy Fe K emission and Fig. 2(b) shows a similar plot for the 
higher energy emission. Consistent with the spectral fitting 
results, it can be seen that the lower energy Fe K emission 
peak is marginally redshifted from the rest-frame energy of 
Fe I Ka (dashed line) at ~ 90% confidence. If we identify 
the higher energy feature with Fe XXVI Lya emission, this 
line is also marginally (at > 90% confidence) redshifted with 
respect to systemic velocity (dashed line). Fig. 3(a) shows 
a two-parameter joint confidence contour of the equivalent 
width (EW) versus FWHM for the lower energy Fe K emis- 
sion and Fig. 3(b) shows a similar plot for the higher energy 
emission peak. 

Fig. 4(a) shows a close-up of the best-fitting Gaussian 
model fit (see Table 1) to the IC 4329 A HEG spectrum. 
Again, the offset of the Gaussian peak energies from the 
dashed lines in Fig. 4(a), which mark the energies of Fe I Ka 
and Fe XXVI Lya emission in the observed frame (6.299keV 
and 6.856keV respectively), show that both the lower and 
higher energy peaks (if Fe I Ka and Fe XXVI Lya respec- 
tively) are indeed marginally redshifted in IC 4329A. If the 
higher energy feature is Fe XXVI Lya emission, then un- 
der some circumstances we should also expect significant 
Fe XXVI Ly/3 emission (e.g. see Bautista et al. 1998; Bautista 
& Titarchuk 1999). We therefore added a third Gaussian com- 
ponent to the model in Table 1 and we extended the HEG 
data energy range to 3.0-8.3 keV since Fe XXVI Ly/3 oc- 
curs at 8.251 keV in the rest-frame (for all successive mod- 



Table 1. Double Gaussian Model 



Fig. 4.— IC 4329A HEG photon spectra between 5.0 and 7.3 keV (ob- 
served energy), binned at 0.01A. The solid lines correspond to physically 
different (and degenerate) models (see text for details). Note that the data are 
not unfolded: the spectra were made by multiplying the ratio of the counts 
in each spectral bin to the predicted model counts in that bin, by the best- 
fitting model. In each case the continuum is a sim ple p ower law and the 
emission-line models are ( a) a double Gaussian (see iQ.ll and Table 1); ( b) a 
single disk line (see i|3.2l and Table 2), and (c) a double disk line (see i!3.4l 
and Table 3). Dashed lines indicate the observed-frame energies expected for 
Fe I Ka (6.400 keV) and Fe XXVI Lya (6.966 keV) emission respectively. 



els we shall return to fitting the data range 3.0-8.0 keV). 
Fe XXVI Ly/3 is hardly detected (AC improves by only ~ 1 
for the addition of a Gaussian with three free parameters, and 
EW< 49eV). The value of the Fe Ly/3:Lya ratio is a power- 
ful diagnostic since it can determine the process of excitation. 
However, the measurement errors on the Fe Ly/3:Lya ratio in 
these data are too large to obtain any meaningful constraints. 
On the other hand, since IC 4329A is one of the brightest 
Seyfert galaxies (sometimes it is brighter than 3C 273), it is 
one of the best candidates to constrain the Fe Ly/3:Lya ra- 
tio with future missions such as Astro-E2. We obtained the 
observed-frame 2-10 keV flux (for comparison with values 
in the literature) for the Chandra HEG spectrum by extrapo- 
lating the above best-fitting model down to 2 keV and up to 
10 keV, and this yielded a value of 16.2 x 10~ n erg cm -2 s -1 , 
corresponding to an intrinsic, source-frame 2-10 keV lumi- 
nosity of 9 x 10 43 erg s" 1 . 

Since we detected an emission line peak that could be due 
to Fe XXVI Lya we searched for evidence of other features 
due to highly ionized Fe (as detected e.g. in NGC 7314 by 
Yaqoob et al. 2003a) by adding a third Gaussian model com- 
ponent. We find EW< 17eV forFe XXV (f) (at the same offset 
velocity and the width of the higher energy Fe K feature in Ta- 
ble 1 ) and we find marginally significant, narrow Fe XXV (f) 
emission at the galaxy systemic velocity (detected at > 95% 
confidence, since AC improves by ~ 5.6 for the addition of 
one free parameter). 
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NOTE. — Model parameters for the best-fitting double-Gaussian model for the 
double-peaked Fe K complex in the HEG spectrum IC 4329A (see i|3.1l . Errors are 
90% confidence for one interesting parameter (AC = 2.706). All parameters are in 
the source rest frame (z = 0.016054). a We show redshifi of the line centroid energy 
as a velocity offset when the low and high energy peaks are identified with Fe 1 Ka 
emission and Fe XXVI Lya respectively. Fe I Ka emission and Fe XXVI Lya emission 
have rest energies of 6.400keV and 6.966keV respectively. Velocities are rounded to the 
nearest 5 km s _1 . Compare these values with the systematic uncertainty in the HETGS 
energy scale of ± 420, ± 460 km s~' at the observed energies of the emission peaks. b 
Gaussian emission-line intensity in units of 10~ 5 photons cm~ 2 s~ L . c Change in the C— 
statistic upon addition of the particular Gaussian model component to a simple power- 
law model for the continuum (there are three additional degrees of freedom for each 
Gaussian component). d Degrees of freedom. 



Next, we replaced the double-Gaussian emission-line 
model described above with a single, broad emission line 
model from a relativistic disk around a Schwarzschild black 
hole (e.g. Fabian et al. 1989). These data cann ot di stinguish 
between Kerr and Schwarzschild metrics (see ^3.6t . The in- 
ner radius of the disk was fixed at 6r g since it could not be 
constrained. The outer disk radius and line intensity were 
free parameters. The radial line emissivity per unit area was a 
power law, r q , where q is the emissivity index. This time we 
also included another disk-line component to model Fe K/3 
emission. Although it makes little difference to the results, 
we included it for completeness. The rest energy was fixed 
at 7.056 keV (corresponding to Fe I emission) for simplicity 
even though the energy of Fe Ka was allowed to vary (but the 
proximity of the latter energy to the Fe I Ka energy, combined 
with the relative weakness of K/3 easily justifies this assump- 
tion). The Fe I K/3 disk-line parameters were coupled to those 
of the lower energy Ka line, with the intensity of K/3 tied to 
(17/150) times the Ka line intensity. There were a total of 
five free parameters for this emission-line model, including 
the intrinsic rest energy of the monochromatic emission line. 
The best-fitting results for this model are given in Table 2 and 
the data and model are shown in Fig. 4(b). The complex Fe-K 
emission can clearly be fitted with a single broad line given: 
a flat radial emissivity (q > -0.72) with emission from ~6- 



70 r g and 



24^j degrees (90%, one-parameter uncertain- 



3.2. Single Relativistic Disk Line Model 



ties). The best-fitting rest-frame energy for this relativistically 
broadened line is 6.74^[J 22 keV, which is indicative of He-like 
Fe. The C-statistic for the single disk-line model is higher (by 
~ 8.5, corresponding to > 99% confidence) than that for the 
double-Gaussian fit in Table 1 . However, the single disk-line 
model we use may be too simplistic (in particular with respect 
to the assumed simple emissivity law and the assumption of 
axisymmetric emission), so we cannot rule it out on the basis 
of the fit statistic alone. Note also that the flat line-emissivity 
index may be a result of fixing the inner disk radius at 6r g . 
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Table 2. Single Disk Line Model and He-like Fe Absorption 
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NOTE. — Model parameters for the best-fitting relativistic Schwarzschild 
disk line models (with and without Fe XXV He-like resonance ab sorp tion) 
fitted to the Fe K complex in the HEG spectrum of IC 4329A (see All 
parameters are in the source rest frame (z = 0.016054). Errors are 90% con- 
fidence for one interesting parameter (AC = 2.706). The inner radius of the 
disk-line model component was fixed at 6r g . " Model fitted with no Fe XXV 
He-like resonance absorption. * Model includes an inverted Gaussian with 
center energy fixed at 6.7 keV to model Fe XXV He-like resonance absorp- 
tion (best-fitting absorption-line parameters are given in the last column of the 
table). c There are local minima all the way up to the rest energy of Fe XXVI 
Lya (6.966 keV). * Fe-K emission-line intensity in units of 10~ 5 photons 
cm~ 2 s~' . e Change in the C-statistic upon the addition of the disk-line model 
component to a simple power-law model, f Degrees of freedom. 



If the line emission is truncated at a larger radius, a steeper 
emissivity law is allowed (e.g. see Done et al. 2000). When 
we fixed q at -2.5 (a value typically found in relativistic Fe K 
line fits to CCD data) and allowed the inner radius to vary, a 
fit was obtained which was only marginally worse (| AC\ < 3), 
and yielded an inner radius of 28! J^. 

3.3. Single Disk Line with Fe XXV He-like Absorption 

It is conceivable that between the low-energy and high- 
energy peaks of the Fe K emission-line complex, we may have 
Fe XXV He-like resonance absorption at 6.7 keV. Such reso- 
nance absorption cutting away the Fe K emission-line pro- 
file has been suggested for NGC 3516 (Nandra et al. 1999) 
and MCG -6-30-15 (Fabian et al. 2002) and discussed in 
a theoretical context by, for example, Ruszkowski & Fabian 
(2002). Therefore, we tested the data against a simple model 
of Fe XXV He-like resonance absorption by adding an inverted 
Ga ussia n component to the single disk-line model described 
in 33.21 The center energy was fixed at 6.7 keV in the rest 
frame of the source, which meant that two more free parame- 
ters were added to the model (the width and equivalent width 
of the absorption line). The results for both the relativistic Fe 
K emission line and the absorption line are shown in Table 2 
so that they can be directly compare d wi th the fit with the rel- 
ativistic Fe K emission line only ( 33.21 It can be seen that 
the addition of the absorption line has little effect on the disk 
inclination angle, the disk outer radius, and the radial emissiv- 
ity derived from the broad Fe K line. Naturally, the intensity 
and equivalent width of the Fe K emission line are required 
to be larger due to the presence of the absorption line. The 
derived equivalent width of the absorption line is 46!^| eV, 

and its width is 12, 1 20^53™ km s" 1 . These values depend 
somewhat on how the underlying emission line is modeled. 
However, since the addition of the absorption line decreases 



the C-statistic by only 6.2 for the addition of two free param- 
eters (corresponding to a significance of < 99%), we do not 
pursue the matter further. 

3.4. Double Disk Line Model 

Next, we fit the Fe K emission-line complex with two Ka 
disk-line model components. The idea here is to investigate 
in a very simplistic way, the possibility of line emission from 
a disk with a vertically stratified ionization structure. For ex- 
ample, the higher energy peak could come from a highly ion- 
ized top layer, whilst the lower energy line emission could 
come from a cooler layer beneath the hot layer. Pairs of pa- 
rameters corresponding to the inclination angle, the inner and 
outer radii for the two disk-line components were tied so that 
they were forced to vary together. Both disk-line model com- 
ponents had inner radii fixed at 6r g . The rest-energy for the 
lower energy line component was a free parameter but the 
rest energy of the higher energy line was fixed at the rest- 
energy of Fe XXVI Lya (6.966 keV). Thus, the radial line 
emissivities of the two disk-line components were indepen- 
dent. The above simplifying assumptions are necessary, given 
the limitations of the data, in order not to over-parameterize 
the model. Again, for completeness we added an extra disk- 
line component to model Fe K[3 emis sion associated with the 
lower energy peak, as described in 33.21 for the single disk- 
line model. However, inclusion of the Fe K/3 emission does 
not add any more free parameters to the model. Thus, there 
were a total of six free parameters in this model of the Fe K 
emission complex. The best-fitting results for this model are 
given in Table 3 and the data and model are shown in Fig 4(c). 
Clearly both disk lines have very similar radial emissivity pro- 
files and are strongly constrained to a near face-on inclination 
(6 < 6.0° at > 95% confidence). The EWs of the low and 
high energy features, as modeled by disk lines, are larger than 
the corresponding EWs in Table 1 (double-Gaussian model) 
and their combined EW is larger than the single disk-line fit 
to the Fe K region (at 90% confidence). The fit statistic is an 
improvement over the double-Gaussian model fit at > 95% 
confidence and is an improvement over the single disk-line 
model at > 3er confidence (AC = -12.6 for one additional free 
parameter). 

3.5. Relativistic Disk Line and Gaussian Model 

Next, we modeled the two peaks in the Fe K band with 
a relativistic disk line for the lower energy emission compo- 
nent and a Gaussian for the higher energy comp onen t. Again, 
we included Fe K/3 emission as described in 33.21 for com- 
pleteness. The best-fitting results for this model are given 
in Table 4. The disk inclination angle is tightly constrained 
to be 12 ±2 degrees and q is flat. The Gaussian best-fitting 
parameters to the high-energy feature are similar to those in 
Table 1, as we should expect. This model is not significantly 
better than the double-Gaussian model statistically (Table 1), 
but it is better than the single broad disk-line fit (Table 2) at 
> 95% confidence since the C-statistic decreases by ~ 10.3 
for three additional free parameters (relative to the single disk- 
line model). Note that the radial emissivity power-law index 
(q) for the disk-line emission is very flat. 

3.6. Complex Continuum 

So far, we have not included the effects of a Compton- 
reflection continuum, or indeed any other complexity in the 
continuum. First of all, we note that there is little evidence 
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Table 3. Double Disk-Line Model 



Table 4. Disk Line Plus Gaussian Line Models 



Low-Ionization Fe K Fe XXVI Lya 



E(keV) 


f. 443+0 050 
D -0.016 


6.966 fixed 


<f 


-2 44+0-22 


-2 20+ 033 


9{°) b 


<6.0 




r out/ r g 


633+ 57 
° -369 




1° 

EW(eV) 


20.6+1* 
120+" 


14 S+ 6 ' 3 

• -7.0 

101+| 


r 


i 7 n+oT8§ 

1 ' /y -0.05 




AC 1 


-25.2 


-22.9 


C (d.o.f.) e 


514.1 (507) 





NOTE. — Parameters for the best-fitting double relativistic disk-line fit to the Fe K 
line region of the HEG spectrum in IC 4329A. The model includes an Fe K /3 emission 
component for the lower-energy line (see i |3.4l for fitting details). Errors are 90% confi- 
dence for one interesting parameter (AC = 2.706) except where noted. The inner radius 
of all disk-line model components was fixed at 6r g . All parameters are in the source rest 
frame (z = 0.016054). The energy of the Fe XXVI Lya disk line component was fixed at 
its rest energy of 6.966 keV. a Power-law index for the line radial emissivity. The lower 
energy Fe Ka and Fe XXVI Lya disk line components were coupled in all parameters 
except this one (q). b Both 9 and r out /r g for the Fe XXVI Lya disk line were forced 
to their corresponding values for the lower energy disk line during the fitting. The C- 
statistic varies non-monotonically close to the best fit. For the disk inclination angle the 
95% confidence upper limit is shown. Formally, 9 < 0.35° at 90% confidence but the 
95% confidence value is more reliable since it is in the region where C is well-behaved. 
F° r r aM jr s , the errors represent 90% confidence values, but only when AC never drops 
below 2.706 again. ' Line intensity in units of 10~ 5 photons cnT 2 s~'. Change in C— 
statistic upon the addition of the particular disk line model component to a power-law 
continuum only. e Degrees of freedom. 



of an Fe K edge above 7 keV in the data (Fig. 4). An Fe K 
edge could be produced either in transmission or reflection 
from optically-thick matter. We adde d a s imple edge model 
to the disk-line plus Gaussian model ( 33. 5> . at 7. 1 1 keV (cor- 
responding to neutral Fe) and found a best-fit value of zero, 
and a 90%, one-parameter upper limit of tf s < 0.054. This 
is comparable with Tp e ~ 0.03 from the XMM-Newton obser- 
vation of IC 4329A by Gondoin et al. (2001). This implies a 
total neutral absorbing column of Ng < 4.9 x 10 22 cm" 2 in the 
line-of-sight. This is certainly consistent with the low-energy 
MEG data (McKernan, Yaqoob, & Reynolds 2004, in prepa- 
ration). 

Due to the restricted bandpass and limited signal-to-noise, 
the data are not sensitive to the Compton-reflection contin- 
uum, which becomes important above ~ 7 keV. However, 
we demonstrated this directly by fitting a model consisting 
of a power law with a high-energy exponential cut-off, and 
Compton reflection from optically-thick, neutral matter (us- 
ing the XSPEC model pexrav, as described in Magdziarz 
& Zdziarski 1995), a Fe Ka line to model the lower energy 
peak in the Fe complex, and another emission line to model 
the higher energy peak. A Fe Kj3 line was also included , for 
completeness, and tied to the Fe Ka line as described in 33.21 
The Fe Ka line, the Fe Kj3 line, and the Compton reflection 
continuum were all blurred with a relativistic kernel based on 
the Kerr metric (Laor 1991; see also Fabian et al. 2002). The 
reason for now using the Kerr metric (for a maximally spin- 
ning black hole), instead of the Schwarzschild metric, is to di- 
rectly demonstrate that the data are not sensitive to the black- 
hole angular momentum. Since the higher energy peak of the 
Fe complex (possibly due to Fe XXVI Lya) in the IC 4329A 
HEG data is quite narrow, it was modeled with a Gaussian and 
not subjected to the relativistic blurring. Also, we included no 
Compton-reflection continuum associated with the the higher 
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NOTE. — Model parameters for the best-fitting relativistic disk line plus Gaussian 
models to the Fe K complex in the HEG spectrum of IC 4329A (see !3.5l for fitting 
details). The model also includes Fe K /3 emission for the lower energy line (see 13.21 
for details). Errors are 90% confidence for one interesting parameter (AC = 2.706). The 
inner radius of the lower energy Fe Ka disk-line model was fixed at 6r K . All quantities 
are quoted in the source rest frame (z = 0.016054). " Simple model with a power-law 
continuum, and a relativistic disk line for a Schwarszchild black-hole. Relativistic 
disk line from maximally rotating Kerr black-hole, with a power-law plus Compton- 
reflection continuum blurred by the relativistic line kernel. Fe XXVI Lya has a rest 
energy of 6.966 keV. L Velocities are rounded to nearest 5 km s" . Line intensity in 
units of 10~ 5 photons cra" 2 s+. e Effective 'reflection fraction' (i.e. R = 1 corresponds 
to the steady-state Compton reflection amplitude from an infinite disk subtending 27T 
solid angle at the X-ray source). ' Degrees of freedom. ? Change in the C-statistic upon 
addition of only the disk line to the continuum. /! Change in the C-statistic upon addition 
of a Gaussian emission-line component to the disk-line plus continuum model. 



energy peak. This is because if it is due to highly ionized Fe, 
the curvature in the reflected continuum is much less than the 
case for a disk which is not ionized, especially for data in the 
restricted bandpass (3-8 keV) that we are dealing with here. If 
the higher energy feature is not a separate line, then the single 
Compton-reflection continuum already included is adequate. 

The free parameters of the laor model were the disk in- 
clination angle, 8, the power-law emissivity index, q, and the 
outer disk radius, r out . Since even the lower energy emission 
peak is not extremely broad, we found the fit to be quite in- 
sensitive to the inner disk radius, so again fixed it at 6r g . The 
Compton-reflection model (pexrav) does not allow the disk 
inclination angle to be tied to that of the laor model, and also 
does not allow inclination angles < 18°. Therefore we fixed 
the inclination angle at 30° for the Compton-reflection con- 
tinuum. This is adequate for our purpose. We also assumed 
solar element abundances and an exponential cut-off energy 
of 200 keV for the intrinsic power-law continuum. This is 
consistent with broadband observations of IC 4329A as well 
as results for other AGNs (e.g. Madejski et al. 1995; Malizia 
et al. 2003; Perola et al. 2002). In any case, our conclusions 
are not sensitive to this assumption. The Compton-reflection 
continuum then had only one free parameter in the fit, the ra- 
tio, R, of the amplitude of the reflected continuum to that ex- 
pected from a centrally illuminated infinite disk, subtending a 
solid angle of 2ir at the X-ray source. In all, this model had 
eleven free parameters, including the overall continuum nor- 
malization. This is one more free parameter than the equiva- 
lent model that utilized the Schwarzschild metric and that did 
not include a Com pton-reflection continuum (i.e. the model 
described in 33.5ll. 
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The results are shown in Table 4 so that they can be directly 
compared with the results when Compton reflection was not 
included, and the continuum was modeled only with a sim- 
ple power law. It can be seen that the quality of the fits for 
the two models, the best-fitting values for the emission-line 
parameters and their statistical error ranges are virtually in- 
distinguishable. This applies to both the relativistic and non- 
relativistic emission lines. It is important to note, however, 
that inclusion of the Compton-reflection continuum does al- 
low the line radial emissivity index to be steeper, relieving 
the requirement that the emissivity law be flat (which may 
be physically difficult to achieve). Another noticeable differ- 
ence between the simple and complex continuum fits is that 
the latter fit yields a somewhat steeper intrinsic X-ray con- 
tinuum power-law index. For the fit with a complex con- 
tinuum, the reflection fraction, R, was consistent with zero, 
and we obtained an upper limit of R < 2.9 (90% confidence, 
one parameter). This confirms that the data are insensitive 
to the Compton-reflection continuum. The result also justi- 
fies using only a simple po wer-l aw for the continuum in the 
model-fitting described in 33.11 - 33.51 Since the relativis- 
tically smeared Compton-reflection continuum is smoother 
and produces less curvature in the net spectrum than an un- 
smeared continuum (in the energy band relevant here), we 
tested the data against the same model with the smearing 
of the Compton-reflection component turned off. Physically, 
such a component could originate in optically-thick distant 
matter, such as a parsec-scale torus. Again, no reflection was 
formally required by the data, with a 90% confidence upper 
limit of R < 1.7. These results further justify using only a 
simple po wer la w for the continuum in the model-fitting de- 
scribed in jni- 

4. DISCUSSION 

In a Chandra HETGS observation of IC 4329 A, we have 
measured broad and complex line emission in the Fe K band 
using the highest spectral resolution currently available in the 
Fe K band (~ 1860 km s" 1 FWHM at 6.4 keV). The line pro- 
file is double-peaked, the peaks occurring at ~ 6.3 keV and 
~ 6.9 keV in the source rest frame. Modeled with Gaus- 
sians, the lower and higher energy peaks have widths of 
20,825!7°^5° km s" 1 FWHM and 3680+g^ km s" 1 FWHM 
respectively. The lower energy peak is resolved by the high 
energy grating (HEG) at > 99% confidence, whilst the higher 
energy peak is only resolved at < 90% confidence. The low- 
energy peak is redshifted from the expected rest-frame energy 
of Fe I Ka by at least 650 km s" 1 and the high-energy peak 
is redshifted from the expected rest-frame energy of Fe XXVI 
Lya by at least 950 km s" 1 . Both of these lower limits on the 
redshift correspond to 90% confidence levels for one parame- 
ter. The emission-line peak at ~ 6.9 keV that we found in the 
Chandra HETGS data has not been detected in any previous 
observations of IC 4329A. We also detected a much weaker 
emission peak at the expected energy of the Fe XXV forbid- 
den line energy at the systemic velocity of IC 4329 A, but the 
feature is detected only with a marginal significance (> 95% 
confidence). 

The amount of Fe in the line-of-sight may be insufficient to 
produce all of the observed EW of the lower energy ^6.3 keV 
Fe K feature (HO^q eV from Table 1) since the measured 
(HEG) upper limit on the optical depth of an Fe K edge cor- 
responds to an absorbing column density < 4.9 x 10 22 cm" 2 
and this is consistent with the low-energy MEG data from the 
same observation (McKernan, Yaqoob & Reynolds 2004 in 



preparation; see also Gondoin et al. 2001). Such a column 
density can at most produce an EW of < 50 e V, even if it fully 
covered the sky (e.g. Krolik & Kallman, 1987; Yaqoob et al. 
2001). To investigate the contribution to the Fe K line emis- 
sion from distant matter, we added an extra (narrow) Gaus- 
sian line component to the two Gaussian model described in 
Table 1. The energy of the new Gaussian component was 
fixed at 6.4 keV in the rest-frame of the AGN since we are 
testing the data for a line which originates from matter suffi- 
ciently distant that there is no gravitational redshift. Also, a 
line originating from such distant matter would be unresolved 
by Chandra , so we fixed the width of the Gaussian at a value 
(1 eV) much less than the instrument resolution. The two 
other Gaussian model components were allowed to vary as 
before. The best-fitting model did not result in a significant 
improvement in the fit-statistic (AC = 2.3 for one additional 
free parameter). Thus, the values and error ranges of all the 
model parameters did not change significantly compared to 
the values shown in Table 1, which were obtained without the 
distant-matter line. The 90% confidence (one-parameter) up- 
per limit on the distant-matter line EW we obtained is 18 eV, 
so that it could contribute at most ~ 1/4 of the low energy 
(6.2-6.4 keV) Fe K fine emission in IC 4329A. This and the 
fact that the low-energy peak is at 6.3 keV and resolved by 
the HEG suggests that the contribution from distant matter is 
not as important as in some other AGN in which the Fe K 
line emission peaks at 6.4 keV and is unresolved by the HEG 
(see Yaqoob & Padmanabhan 2004). Thus, it is likely that 
the bulk of the ~ 6.3 keV line emission originates in an ac- 
cretion disk. It is difficult to deduce what the ionization state 
of the matter must be because of the possibility of a peak at 
energies higher than 6.4 keV being gravitationally redshifted 
down to ~ 6.3 keV. Note that a radius of ~ 65r g is required to 
gravitationally redshift a 6.4 keV line to 6.3 keV. 

The interpretation of the overall, double-peaked Fe K line 
profile is difficult. Several physically distinct models can de- 
scribe the data. The two peaks may correspond to the Doppler 
horns of a single disk line or they may be artificial in the sense 
that He-like Fe K absorption has removed the central portion 
of a disk line profile which may not have prominent Doppler 
peaks. Alternatively, the higher energy peak at 6.9 keV may 
correspond to Fe XXVI Lya emission which is independent of 
the line emission associated with the lower energy peak. In 
that case, the Fe XXVI Lya line may or may not originate in a 
disk. Again the two scenarios are degenerate as far as the data 
are concerned. If there is Fe XXVI Lya line emission from a 
disk it may come from a layer which is hotter than that pro- 
ducing the lower-energy line (for example the hot corona, e.g. 
see Nayakshin & Kallman 2001; Ballantyne, Ross, & Fabian 
2001). Although the ~ 6.9 keV line appears to be fairly nar- 
row we have shown by direct model fitting that it could still 
originate from a disk which extends down to 6r g , with a radial 
line emissivity per unit area which is not flat. 

The line emission peak at 6.9 keV is reminiscent of a sim- 
ilar feature found in a long XMM-Newton exposure of the 
Seyfert galaxy MCG -6-30-15 (e.g. Fabian et al. 2002). 
In that case there is also an ambiguity between Fe XXV He- 
like absorption and and Fe XXVI Lya emission, and the peak 
energy is slightly redshifted from the expected energy of 
Fe XXVI Lya, just as in IC 4329 A. Fe XXV He-like absorption 
has actually been observed with less ambiguity in some AGN 
(PG 1211+143, Pounds et al. 2003a; PG 0844+349, Pounds 
et al. 2003b; NGC 3783, Reeves et al. 2003). 

Although there are several interpretations of the complex 
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Fe K emission revealed by the HEG spectrum of IC 4329 A, 
each of the different scenarios is interesting from the point of 
view that future observations with Astro-E2 will break much 
of the degeneracy and whichever scenario turns out to be cor- 
rect, we would learn something new about the central engine. 
In particular, observational constraints on the ionization struc- 
ture of the accretion disk in AGN are currently very poor and 
IC 4329A may be a good astrophysical laboratory for improv- 
ing them. 
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